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We revisit the problem of constraining steps in the inflationary potential with cosmological data. 
We argue that a step in the inflationary potential produces qualitatively similar oscillations in the 
primordial power spectrum, independently of the details of the inflationary model. We propose 
a phenomenological description of these oscillations and constrain these features using a selection 
of cosmological data including the baryonic peak data from the correlation function of luminous 
red galaxies in the Sloan Digital Sky Survey. Our results show that degeneracies of the oscillation 
with standard cosmological parameters are virtually non-existent. The inclusion of new data severely 
tightens the constraints on the parameter space of oscillation parameters with respect to older work. 
This confirms that extensions to the simplest inflationary models can be successfully constrained 
using cosmological data. 
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I. INTRODUCTION 

Recent data from the Wilkinson Microwave Anisotropy 
Probe (WMAP) % i, i, i] observations of the 
anisotropics of the Cosmic Microwave Background 
(CMB) are in excellent agreement with the predictions 
of inflationary cosmology. In its simplest implementa- 
tion, inflation is driven by the potential energy of a single 
scalar field slowly rolling down the potential towards the 
real vacuum. Under the assumption that the potential 
is sufficiently flat and smooth, the resulting spectrum of 
density perturbations is almost scale-invariant and can 
be described with a power-law. In the context of this 
slow-roll paradigm, a number of authors have used the 
WMAP data and various other complementary data sets 
to derive bounds on the inflationary parameter space. 
These include constraints on specific inflationary models 
H 00,11, the Hubble dynamics during inflation [j| [l(| , 
or, in a more empirical fashion, the parameters charac- 
terising the primordial power spectrum [111 . Il2t Il3l . |14| . 

In more general classes of inflationary models, however, 
slow roll may be violated for a brief instant Jig, [l6[ • In 
single-field inflation models, such an effect can be mod- 
elled by introducing a feature such as a kink, bump or 
step [17| to the inflaton potential. A step, in particular, 
can be regarded as an effective field theory description 
of a phase transition in more realistic multi-field models 
fl8l |. which may arise naturally in, e.g., supergravity- [l9[ 
or M-theory-inspired inflation models [2fj |. 

This interruption of slow-roll will leave possibly de- 
tectable traces in the primordial power spectrum. Specif- 
ically, wavelengths crossing the horizon during this fast- 
roll phase will be affected [2lT [22J, leading to a devia- 
tion from the usual power-law behaviour at these scales. 
Such non-standard power spectra have been brought for- 
ward to explain the peculiar glitches in the temperature 



anisotropies 12311 a s well as the observed low power at the 
largest scales [24. l25j . 

Step-like features in the inflaton potential will lead to a 
burst of oscillations in the primordial power spectrum. A 
particular realisation of a step potential was confronted 
with the data in Ref. [26fl for fixed cosmological parame- 
ters and more generally in Ref. (27| . 

In the present work, we extend the analysis of [27| 
in several important aspects. Firstly, we generalise our 
method to spectra corresponding to a whole class of step- 
inflation models with arbitrary (slow-roll) background in- 
flaton potentials. This allows us to derive constraints on 
parameters characterising the feature in a more model- 
independent way. Secondly, we address the question of 
parameter degeneracies: could the presence of a feature 
bias our estimates for the values and errors of the cos- 
mological parameters, such as the baryon or dark matter 
density, in any way? Thirdly, we consider new data- 
sets: apart from CMB and matter power spectrum data 
sets, we consider also the measurements of the position of 
acoustic peak in the real space two-point galaxy correla- 
tion function data (BAO) from the luminous red galaxy 
(LRG) sample of the Sloan Digital Sky Survey (SDSS) 
[28|. This data is, in principle, especially well suited to 
constraining (or detecting) small amplitude oscillations 
in the power spectrum which would show up as a peak 
in the correlation function. However, due to biasing and 
weakly non-linear structure formation, this data is dif- 
ficult to interpret and we pay special attention to do it 
carefully. 

The paper is organised as follows: In Sec. [Til we briefly 
remind the reader of the exact formalism of calculating 
the power spectrum from a given inflaton potential and 
compare it with the slow-roll approximation. In Sec. IIII1 
we discuss the dynamics of the inflaton field rolling over a 
step and introduce a generalised step model. Section [IV] 
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is dedicated to our analysis methods with an emphasis 
on the determination of the likelihood for the BAO data 
set. We present our results in Sec. [V] and, finally, draw 
our conclusions in Sec. IVIl 



II. INFLATIONARY PERTURBATIONS 

Let us start this section with a brief recapitulation of 
how to calculate the primordial spectrum of curvature 
perturbations V-r.-, using the formalism of Stewart and 
Lyth HI. 

In the following, we will set c = h = 8irG = 1. We 
consider the gauge invariant Mukhanov variable u [30l.[3l| 
given in terms of the curvature perturbation 1Z: 



-zK. 



(1) 



For our purposes, it is convenient to introduce another 
time parameter, the number of e-foldings, defined by 
N = In a. In terms of N, Equations (|2|), and read 
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(9) 
Uk = 0, 



with Nq determining the normalisation of the scale fac- 
tor. This coupled system of differential equations can 
easily be solved numerically, once a suitable set of initial 
conditions has been chosen. 



Here, z = a<f>/H, where a is the scale factor, cf> the infla- 
ton field, H the Hubble parameter and the dot represents 
a derivative with respect to time t. The Fourier compo- 
nents of u obey the equation 



Uk = 0, 



(2) 



with a prime denoting a derivative with respect to con- 
formal time r. 

Finally, we can define the primordial power spectrum 
of curvature perturbations Viz{k) via the two-point cor- 
relation function 



2ttJ 

k 3 



Vn(k) <5( 3 >(fci-fc 2 ). 



(3) 



Assuming gaussianity and adiabaticity, this quantity con- 
tains all the necessary information for a complete statis- 
tical description of the fluctuations. It is related to Uk 
and z via 



2^2 



z 



(4) 



A. Background Equations of Motion 



In order to find a solution to Equation ©, one needs 
to know the behaviour of the term z" / z. Its evolution is 
determined by the dynamics of the Hubble parameter and 
the unperturbed inflaton field, governed by Friedmann's 
equation 



H 2 



(v + U 2 ) 



and the Klein-Gordon equation for (j) 

■ dV 

4> + 3H(j)+ — = 0. 



(5) 



(G) 



B. Initial Conditions 

Supposing that at a time N sr the system has reached 
the inflationary attractor solution 



and is rolling slowly, 



2 « V ((/>), 



(10) 
(11) 



the initial conditions for <fi and H will be given by 

<t>{N sr ) = <t>sr (12) 

1 dV 



4>,N(N sr ) 



H(<j> sr ) 



v(<p s 



(13) 



(14) 



The initial conditions for Uk can be obtained by requiring 
the late time solution of © to match the solution of a 
field in the Bunch-Davies vacuum of de Sitter space, given 
by 



Uk(r) 



2k 



kr 



(15) 



at early times, well before the observationally relevant 
scales leave the horizon. For k ^> z" / z (or, equivalently, 
kr 3> 1) this can be approximated by the free field solu- 
tion in flat space 

u k = -^=e- ikT . (16) 



Fixing the irrelevant phase, we obtain the initial condi- 
tions for a mode k 



Uk(T Q ) 



1 



2k 



u 'k( T o) = -*'y 2 
at a time tq satisfying k 3> z" /z\ To . 



(17) 
(18) 
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C. Slow Roll 

Following Ref. [Hj], we define the Hubble slow roll 
parameters by 



'V? 



H 



n 



d In H « 
din a 



l/n 



(Jj(l))n-ljy(n+l) 



H n 



l/n 



(19) 

for n > 1, with a superscript "(n)" denoting the nth 
derivative with respect to (j). In addition to that, we 
define °(3h = 2{H^ / H) 2 . The first three parameters of 
the Hubble slow roll hierarchy read 



e H = °(3 H = 2 
r]H = 1 0H 



H{4>) 



H 
H 2 ' 



(j)H 



,2 _,2 fl ^2 9 i/ (1) (^ (3) (0) 



H 2 W 



2H 2 (f> 
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(20) 
(21) 

I- 

(22) 



Using these definitions it can be shown that the mode 
equation ^ can be written as 



u'l + r - 2a z H 2 



26^7^ + |77| + £ff ] I Wfc = 



(23) 



Note that this expression is exact: it does not assume 
the slow roll parameters to be small. 

From a model-building point of view, where one re- 
gards the Lagrangian (or the scalar potential) of the the- 
ory as the fundamental quantity, the calculation of the 
Hubble slow roll parameters can be quite involved. In 
this sense it may be more convenient to work with the 
potential slow roll parameters instead, which use deriva- 
tives of the potential instead of derivatives of the Hubble 
parameter. The first three potential slow roll parameters 
are defined by 



v ' 

V 2 



(24) 

(25) 
(26) 



If the attractor condition (Eq. JTUJ) ) is satisfied, the two 
are approximately related via [33j 



^e 2 



e H = e 



0(3), 
V+3 

Ser] + 3e 2 +0(3), 



€ H = e 

Vh = V - e + |e 2 + \eri , 3 

t 2 Q I o,2 



•0(3), 



(27) 
(28) 
(29) 



up to corrections of third and higher orders. Expressed 
in terms of the potential slow roll parameters, z" / z is 
given by 



Z e 2_ 



=2a 2 H 2 [l + fe- 

fer 1+ ^ 2 + ^ 2 +0(3)} 



(30) 



It is commonly assumed that the first two slow roll 
parameters vary slowly with time (i.e., (^h) ^ ■"•)■ Then 
it follows that, if one wants to sustain inflation for long 
enough to solve the horizon and flatness problems, £(#) 
and |?7(_fr)| will also have to be much smaller than unity. 
In this ("slow roll") limit, we have z" / z w 2a 2 H 2 , H « 
and a oc exp[_fft]. 

Let us now turn back to Eq. ^j, which is basically 
the equation of an oscillator with a time dependent mass 
term, and discuss its solutions. The initial conditions 
imply that for wavenumbers with 

- » H, (31) 
a 

i.e., with wavelengths much smaller than the horizon, the 
solution is given by Equation (fT5j) and Uk describes a cir- 
cular motion in the complex plane. Due to the exponen- 
tial growth of the scale factor, the physical wavelenghts 
will be blown up and leave the horizon, eventually sat- 
isfying k/a <^ H. In this limit, the growing solution for 
Uk is given by 

Uk oc z. (32) 

Hence, the spectrum V-ji will converge to a constant 
value for super-Hubble modes, i.e., the perturbations 
"freeze in" . We can also conclude that the fate of a per- 
turbation with wavelength k is decided when k/a ~ H 
and the spectrum will have its final shape imprinted on 
horizon exit. It is not until much later, when the modes 
reenter the horizon during radiation or matter domina- 
tion, that they will exhibit dynamical behaviour again. 

Generically, the power spectrum will not be scale inde- 
pendent, with a scale dependence being induced by the 
variation of, e.g., the potential energy and the Hubble 
parameter as the inflaton field rolls down the potential. 
In the slow roll regime, however, the scale dependence is 
rather weak and Vn can be reasonably well approximated 
by a power law: 

ns — 1 

(33) 



k 



with the normalisation A 
As 



s given by 
1 V 



24tt 2 



and the spectral index 



ko—aH 



1 - 6e + 2r). 



(34) 



(35) 



Before we talk about relaxing some of the assumptions 
that went into this analysis, let us quickly mention infla- 
tionary tensor perturbations. 
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D. Tensor Perturbations 

Apart from the scalar perturbations described above 
inflation also generates tensor perturbations, with a spec- 
trum given by 

fc 3 



T^grav (&) 

and the mode equation 



2vr 2 



a 
a 



v'k + [k 2 -—) v k = 0. 



(36) 



(37) 



This equation is very similar to the scalar one. This 
similarity can be readily seen if we express the "mass 
term" a" /a in terms of the slow roll parameters: 



a 
a 



2a 2 H 2 [1 - ±e H ] (38) 

~ 2a 2 H 2 [1 - ±e + |e 2 - \er) + 0(3)] . 

Just like the scalar modes, tensor perturbations will 
also freeze in at horizon exit. In the slow roll case their 
spectrum is approximately 



with the tensor spectral index 

tit — — 2e, 

and normalisation 

2 



3tt 2 



V 



(39) 
(40) 
(41) 



kQ—aH 



III. SLOW ROLL INTERRUPTED 

The validity of the power-law parameterisation of the 
primordial spectra rests on the assumptions that the 
slow-roll parameters are small and change slowly with 
time. Let us relax the latter and allow e and 77 to change 
significantly on a timescale AN < 1. This has the con- 
sequence that we can also allow e and/or r] to become of 
order unity momentarily, provided that at a later time, 
the system returns to the slow roll regime. We also as- 
sume here that the system starts in a state where the slow 
roll conditions are fulfilled, in order to give it enough time 
to reach the inflationary attractor solution. 

This effect can be modelled by adding a local feature, 
such as a step or a bump, to an otherwise flat inflaton 
potential. 



A. Chaotic Inflation Step Model 

Let us examine the consequences of such a feature us- 
ing as an example the same model potential as in [27] ] 



V(^) = \m 2 ^ 2 M +, 



: tanh 



-b 



(42) 
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FIG. 1: Top: z"/z divided by a 2 H 2 for b = U,c= 10" 3 and 
d = 2 x 10 -2 versus the number of e-foldings. TV is set to zero 
for <f> = b. It takes the inflaton field roughly half an e-folding 
to roll over the step. 

Bottom: Hubble slow roll parameters at the step, eh (dotted 
green line) remains negligible throughout, while rju (solid red 
line) and £h (dashed blue line) violate the slow roll conditions. 



This potential describes standard m 2 <j) 2 chaotic inflation 
[33 | with a step centered around (f> = b. The height of 
the step is determined by c, its gradient by d. We do 
not want inflation to be interrupted by the step, so we 
stipulate \c\ <C 1 to ensure that the potential energy will 
always dominate over the kinetic one. 

As pointed out above, the eventual spectrum crucially 
depends on the dynamics of z" /z, which can easily be de- 
duced from the solution of Equations ([7]) and §8§ . For a 
typical choice of parameters, we plot the numerical solu- 
Generically, we find that z" / (za 2 H 2 ) 



tion in Figure l(a 



has a maximum before the inflaton field reaches b, a min- 
imum shortly afterwards and it will return to the asymp- 
totic slow roll value of ~ 2 after 0(1) e-folding. Com- 
parison with the Hubble slow roll parameters (Fig. 1(b) ) 



shows that this behaviour is mainly caused by rju and 
while en remains small. This is a consequence of the 
condition c <C 1. Beware that the potential slow roll 
approximation Eq. (|30p will in general not work for this 
potential since the contribution of higher derivative terms 
can be large. The smallness of en (and hence e) also im- 
plies that there will not be any sizable deviations from a 
power law for the spectrum of tensor perturbations. 
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So, how will this particular behaviour of z" / z influence 
the solution for Uk and eventually the spectrum compared 
to a model with no step? It is obvious that modes with 
k 2 3> Max\z" /z\, i.e., modes that are well within the 
horizon at the time of the step, will not be affected at all 
and Uk will remain in the oscillatory regime. For k 2 < 
Max|z"/z|, the maximum in z"/z will result in a boost of 
exponential growth for Ufe, reverting to oscillations when 
z" / z goes negative and eventually return to the growing 
solution. We depict the motion of Uk in the complex 
plane in Fig. [5] 




■3-2-10123 -3-2-10123 



lm[u k ] lm[u k ] 

3 i 1 1 1 1 1 1 




-3 1 1 1 1 1 1 1 

-3-2-10123 

lm[u k ] 

FIG. 2: These figures show the evolution of ut in the complex 
plane, where Uk has been normalised to one in the oscillating 
limit. The choice of initial conditions (|17ll8[l ensures that 
the motion will be initially circular. The top left plot shows a 
mode that is not affected by the feature, so that the circular 
oscillation goes straight into a growing motion. In the other 
two plots the circle gets deformed by an intermittent phase 
of growth triggered by the peak of — , to be followed by an- 
other phase of elliptic oscillations (caused by the dip of — ) 
until finally the modes leave the horizon and start growing. 
Whether a mode is suppressed or enhanced by this mecha- 
nism depends on the phase of the oscillation when the growth 
sets in. Growth along the semi-major axis will lead to an en- 
hancement (top right), whereas growth along the semi-minor 
axis entails a suppression (bottom) with respect to the modes 
of the corresponding featureless model. 

When an oscillatory phase is preceded by a growing 
phase, the initial circle will be distorted to an ellipse. As 
the growth sets in again, the the mode will be suppressed 
or enhanced, depending on the phase of the oscillation, 
which itself is ^-dependent. In the spectrum, this can be 
observed as oscillations. This mechanism will be most 
effective for modes that are just leaving the horizon, for 
modes with k 2 <§; Max\z" / z\ the phase difference will be 



negligible. 

Hence, a localised feature in the potential will lead to 
a localised "burst" of oscillations in the spectrum (see 
also Ref. [HI), while large and small scales will remain 
unchanged with respect to the spectra of the asymptotic 
background models. This is shown in Fig. [3) Note that 



cc 
Q_ 




0.1 1 10 100 
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FIG. 3: Primordial power spectrum for a model with 
m = 7.5 x 10" 6 , b = 14, c = 10" 3 and d = 2 x 10" 2 (solid 
black line) with wavenumber k given in units of aH\$=b- The 
dotted red line depicts the spectrum of the same model with 
c set to zero. 

the wavelengths affected by the feature are those that are 
about to leave the horizon as the inflaton field reaches the 
centre of the step. In particular, also the frequency of the 
oscillations of the spectrum is proportional to this scale. 

What remains is to identify the horizon size at the 
step with a physical scale today. This connection can be 
made if one knows the total number of e-foldings A* of 
inflation that took place after a known physical scale fc* 
left the horizon. Technically, we evolve the background 
equations (|7l8p until the end of inflation N en d, (defined 
by a{N ea d) — 0). The scale k* can then be determined 
in units of aH\ r f >= } > via 

, a{N cnd - N*) H(N end - A*) 



As long as the spectrum of the c — background model 
is only mildly scale dependent, there will be a strong 
degeneracy between A* and b: shifting the feature in the 
potential will have the same effect as shifting the scale of 
k. In the following we will therefore not treat A* as a free 
parameter, but set A* = 50 for fc* = 0.05 /iMpc -1 . If 
we want the feature to affect scales that are within reach 
of current observations, this will require b to lie in the 
interval 14 < b < 15. 
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B. Model Dependence 

Having analysed a specific example in the previous sub- 
section, let us now address the question of model depen- 
dence: Will we arrive at different conclusions if we modify 
the background inflationary model (e.g., X(f> 4 instead of 
m 2 4> 2 ) or the parameterisation of the step? 



We will argue that a more general potential 
V{4>) = V + /(<£) 5(0 - b) 



(44) 



leads to a qualitatively similar spectrum as the potential 
(gSJ. Here, Vb S (<f>) = V + /(0) is the background po- 
tential, which should fulfil the slow roll conditions with 
/ and Vo positive definite. The function S((f>) parame- 
terises the step, and should monotonically asymptote to 
1 ± c (c <C 1) for <j) ^> b and cf> <C b, respectively, with 
5(0) = 1. 

As we have seen above, the derivatives of the potential 
are crucial to determining the spectrum. In general, the 
derivatives of V are given by 



^ (n) w = £(f)/ (i) w 

i=0 



i) 



(0). (45) 



Far away from the step, the derivatives of 5 will be 
negligible and the potential and its derivatives are ap- 
proximately 



7(0)~F o + /(0)(l±c)^F b g(0), 



(46) 



F("'W^/W(^)(l±c)^FW(# (47) 

Since the slow roll conditions hold here, the spectrum 
will be given by Eq. (|33[) with 



A? 



ns ~ ng S ± c 



2Vq 
V +f 



In the special case Vo 
A s = Ag s (l ±c) and n s = J 



2 J + <D{c 2 ) 



(48) 
(49) 



0, we have exactly 
If Vq 7^ 0, there are addi- 



tional corrections of order c to the normalisation and also 
corrections to the tilt, which are suppressed by c and the 
slow-roll parameters of the background model. In both 
cases, one asymptotically recovers the spectrum of the 
background model in the limit c -C 1 . 

Near the step, however, the derivatives of V will have a 
contribution from the derivatives of 5. If the step is sharp 
enough, the nth derivative of V will be dominated by the 
nth derivative of 5, since the other terms are suppressed 
with factors of the order of the slow roll parameters of 
the background model. Hence, the dynamics of z" / z near 
the step hardly depends on the background, but is deter- 
mined by the form of 5. On the other hand, any 5 that 
gives a z" / z which roughly shows a behaviour like the 
one depicted in Fig. 1(a) will lead to a burst of oscil- 



lations in the power spectrum. The similarities between 



spectra of different background models are illustrated in 
Fig. [3J where we plot the spectra of a hybrid inflation 
type potential 



V(4>) = V + | m 2 2 ( 1 + ctanh l / 



- b 



, (50) 



and another monomial potential with a different form of 
the step function 



V{(j)) = \(f> 4 I 1 + carctan 



(51) 



Note that despite the difference in background models 
and step functions, the maxima and minima of the oscil- 
lations occur at the same wavelengths. 
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FIG. 4: Primordial power spectra of a hybrid inflation type 



step model (|50j with V = 3.7 x 10 , m = 3.2 x 10~ s , 
b = 0.0125, c = 10" 3 and d = 5 x 10~ 5 (dashed blue line), 
and of potential (|51|) with parameters A = 6 x 10~ 14 , b — 21, 
c = 5 x 10" 4 and d = 0.02 (dotted green line). The hybrid 
inflation background model has ns > 1, suppressing large 
scale fluctuations, while the A</> 4 model has ns < 1 with more 
power on large scales. 

To alleviate the model dependence of the analysis when 
confronting theory with experiment, we choose a phe- 
nomenological approach and define the spectrum of a 
generalised step model 



gsm 



?V tep f 



k 



"S-n| tcp 



(52) 



Here, P n stcp is the spectrum obtained from the poten- 
tial Eq. (|42]) and ng tcp = 0.96 is the spectral index of 
the ■^m 2 (f> 2 model. The quantity ns then describes the 
overall effective tilt of the spectrum. Spectra of this type 
will arise from potentials of the form Eq. (|4"4")l . While 
the fine details of particular models may differ slightly 
from this approximation, Eq. (|52[) will nevertheless cap- 
ture the broad features of a large class of background 
models, since, as argued above, the shape of the burst 
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of oscillations is largely independent of the background 
model. Minor differences would likely be washed out in 
the angular power spectrum of the CMB anyway [3^ |. 
The asymptotic behaviour of models with Vq = will be 
reproduced exactly; for Vq > it will be approximate, 
with errors of order c. 

There is a catch however: in this analysis the parame- 
ters b, c and d will be bereaved of their meaning as param- 
eters of the potential. Instead, they should be interpreted 
as phenomenological parameters which describe the spec- 
trum. This does not preclude us from deriving meaning- 
ful constraints, though. We argued that the shape of the 
modulation of the spectrum is largely independent of the 
background, so similar modulations should be the con- 
sequence of similar step dynamics. A useful quantity in 
this context is the maximum value the slow roll parame- 
ters e, r\ and £ 2 can reach at the step. For the potential 
(j42| . we can estimate e max , ry ma x and ^ M in terms of b, 
c and d: 



2d? + bd' 



^ + 0.77-, 



I Sn 



J d 4 



bd 3 ' 



(53) 
(54) 

(55) 



assuming c < 1, d < 1 and b > 1. Note that £ 2 = for 
the background model. 

Along the same lines, one can can replace b with k s , 
corresponding to today's wavenumber of the perturba- 
tions that left the horizon during inflation when 4> — b. 



IV. DATA ANALYSIS 



C. Generalised step model, which uses an effective tilt 
ns in addition to the parameters of model B. Con- 
straints on e max - e bg , r? max - rj bs and £ 2 lax are 
derived using Eqns. (|53"j) - (|55|) . 

We limit our analysis to scalar perturbations. While ten- 
sor perturbations may, in principle, give a subdominant 
contribution, their spectrum will be smooth in the class 
of models studied here, so we do not expect any major 
degeneracies with the step parameters. 



B. Data Sets 

To assess the influence of different data on the con- 
straints, we perform the analysis for each of the models 
using three different sets of data: 

1. WMAP three year temperature and polarisation 
anisotropy datap], S S, 3 (WMAP3). The like- 
lihood is determined using the October 2006 ver- 
sion of the WMAP likelihood code available at the 
LAMBDA website [H]. 

2. WMAP3 plus small scale CMB tempera- 
ture anisotropy data from the ACBAR [39| . 
BOOMERANG [Hi and CBI [U experiments, 
plus the power spectrum data of the luminous red 
galaxy sample from the Sloan Digital Sky Survey 
(SDSS), data release 4 [42j. To avoid a dependence 
of our results on nonlinear modelling, we only use 
the first 13 fc-bands (k/h < 0.09 Mpc^ 1 ). 

3. Same as data set 2, plus two-point correlation func- 
tion data from the SDSS LRG [H. 



We compare the theoretical predictions of three the- 
oretical models (A, B and C) with observational data. 
We use the Markov chain Monte Carlo (MCMC) pack- 
age cosmomc [37| to reconstruct the likelihood function 
in the space of model parameters and infer constraints 
on these parameters. 



A. Models 

The three models have four parameters in common: 
LOb (baryon density), lo c (CDM density), r (optical depth 
to reionisation) and 9 S (sound horizon/ angular diameter 
distance at decoupling). The difference lies in the pri- 
mordial power spectrum. 

A. Vanilla power-law ACDM model: the initial spec- 
trum is parameterised with A§ and ns- 

B. Step model (Eq. (|4"2"| ) with parameters As, b, c and 
d. 



C. Analysis 

Our constraints are derived from eight parallel chains 
generated using the Metropolis algorithm [43[. We use 
the Gelman and Rubin R parameter [44| to keep track 
of convergence of the chains, stopping the chains at 
R — I < 0.05. Since the likelihood function is highly non- 
gaussian in some parameter directions and even multi- 
modal in certain cases, we double-check our results by 
comparing with chains generated with a variation of the 
multicanonical sampling algorithm [45j . 

D. Priors 

Apart from the hard-coded priors of cosmomc on Hq 
(40 km/s/Mpc < H a < 100 km/s/Mpc) and the age 
of the Universe (10 10 a < A < 2 x 10 10 a), we impose 
flat priors on the other cosmological parameters. For 
the parameters of the potential we choose a flat prior 
on b € [14, 15] and logarithmic priors on c, d and c/d 2 
(logce [-6,-1], logde [-2.5,-0.5], log c/d 2 e [-5,3]). 
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E. Baryon Acoustic Peak 

Oscillations in the dark matter power spectrum due 
to acoustic oscillations in the plasma prior to decoupling 
result in an single peak in the two-point correlation func- 
tion of the distribution of galaxies £(r). In Ref. (28[, the 
authors claim the detection of such a peak and identify it 
as corresponding to the baryonic oscillations of the mat- 
ter power spectrum. 

Since any oscillation of the spectrum, regardless of its 
origin, will lead to a feature in the correlation function, 
this data set is particularly well suited to constraining os- 
cillations in the initial power spectrum as well, provided 
that the features are not completely washed out through 
subsequent evolution. 

The correlation function is related to the matter power 
spectrum P(k) via a Fourier transform: 



Jo kr 



(56) 



Technically, the upper limit of the integral would be some 
ultraviolet cutoff kjjy , chosen such that the error in £ is 
small (<C 1%). For the scales covered by the SDSS data, 
i.e., comoving separations between 12 and 175 /i _1 Mpc, 
this requires a momentum cutoff kuv > 1 h/Mpc. At 
these wavenumbers, however, nonlinear effects cannot be 
neglected anymore, which makes the theoretical predic- 
tion of £ somewhat tricky. 

The standard procedure is outlined in section 4.2 of 
Ref. [28| and involves corrections for redshift space dis- 
tortion, nonlinear clustering, scale dependent bias, and a 
smoothing of features on small scales due to mode cou- 
pling. All of these methods were calibrated with nonlin- 
ear simulations in a vanilla cosmology setting and it is 
not obvious that they should be applicable to our case. 
With the exception of the smoothing, however, the effect 
of these corrections on the correlation function is smaller 
than 10% and will only be noticable at scales < 40/i _1 
Mpc (see Figure 5 of |28j|). So even if we assume a large 
uncertainty in the nonlinear corrections, the accuracy of 
the theoretical correlation function will still be of order 
a few per cent, that is smaller than the error bars of the 
data. 

Let us look at the smoothing procedure in a bit more 
detail. In the usual case, the dewiggled transfer function 
Tdw is a weighted interpolation between the linear trans- 
fer function Tn n and the Eisenstein-Hu [46J no-wiggle 
transfer function T nw 

T dw (k) = w(k) T Un (k) + (1 - w(k)) T nw , (57) 

with a weight function w(k) = exp [— (afc) 2 ] and 
a = IhT 1 Mpc. This is related to the dewiggled spec- 
trum by 

P dw (k) = kTi w (k)V lz (k). (58) 

In the case of a non-smooth primordial power spectrum 
V-ji(k), one should of course also dewiggle the initial fea- 
tures. In order to recover the standard procedure for 



power-law spectra, we will instead smooth the quantity 
f(k) = (P(k)/k) 1/2 = T(k)^lMk). (59) 

The use of the no-wiggle transfer function rests on the 
assumption that at small scales, mode coupling will to- 
tally erase all structure, which is reasonable as long as 
the amplitude of features is of the same order as that 
of the baryon oscillations. For much larger oscillations, 
mode coupling might not be efficient enough to erase all 
structure; it is likely that some residual oscillations will 
remain. So instead of a no-wiggle T nw , we will use a 
smoothed T s defined by 



f s (k,q) = exp 



lnfe+g/2 



dlnfc' In 



1 Jlnk-q/2 



fun(k') 



(60) 



i.e., a convolution of T\ ln with a top hat function of width 
q in log- log space. The dewiggled power spectrum is then 
given by 

Pdw(fc, q) = k (w(k) Tii„(fc) + (1 - w(k)) f s {k, q) 

' (61) 

Without turning to TV-body simulations it would be 
hard to estimate how much the spectrum will have to 
be smoothed, though. Therefore, we will determine the 
BAO likelihood £bao by marginalising over q: 



£bao = J dq£(q)n(q). 



(62) 



V. RESULTS 

An important question in the context of a model- 
dependent analysis is how the choice of model will affect 
the estimates of the parameters, particularly if the mod- 
els are nested. Possible degeneracies between "standard" 
and newly introduced parameters can bias means as well 
as errors. In Fig. [5l we plot the marginalised likelihood 
distributions for the vanilla parameters for all three mod- 
els with data set 1. There are small differences between 
models A and B for flbh 2 , r and the normalisation. These 
arise due to the fact that in model B, the tilt of the spec- 
trum is fixed. There is a well-known degeneracy between 
these parameters and the spectral index. Fixing the tilt 
near the best fit value will reduce the errors on the pa- 
rameters it is degenerate with, which is precisely what is 
happening here. 

The distributions for models A and C show a remark- 
able similarity which leads us to conclude that the pres- 
ence of a feature will not have any statistically significant 
influence on the results for the parameters of the vanilla 
model. This conclusion remains unchanged if we consider 
the other data sets. 

Another interesting question is whether the data prefer 
the presence of a feature over a smooth spectrum. How 
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FIG. 5: Marginalised likelihoods for model A (solid red line), 
model B (dot-dashed line) and model C (solid black line). 
The differences between the results for A and C are marginal. 
For some parameters, the results for model B differ slightly. 
This should be attributed to the degeneracies of the spectral 
index with these parameters and the fact that the tilt of the 
spectrum is fixed in this model, it is not due to the presence 
of a feature. 



much will a feature improve the fit and can we understand 
why? 

In Ref . , we studied model B and found two regions 
in parameter space which improve the fit to the WMAP3 
data by A\ 2 ~ 5 and Ax 2 ~ 7, respectively. The former 
corresponds to oscillations at large scales (I ~ 20 — 30) , 
while the latter has oscillations of a wavelength similar to 
the baryonic acoustic oscillations and lies near the third 
peak of the CMB temperature power spectrum. Adding 
small scale CMB data and, in particular, the power spec- 
trum data of the 2003 data release of the SDSS [47| , im- 
proved the fit of the small scale maximum to A\ 2 ~ 15. 

In the present work, we replaced the old main sample 
data with the luminous red galaxy sample of the most re- 
cent SDSS data release. With this newer data set, how- 
ever, we do not find such an enhancement of the A^ 2 
anymore. In fact, it appears to disfavour a large feature 
near the third peak. Given the better quality of the LRG 
power spectrum data and the fact that the BAO data 
also does not seem to support this effect, it is likely that 
the improvement in the fit was just a fluke. The dis- 
appearance of this maximum of the likelihood function 
is illustrated in Fig. [51 where we show the mean likeli- 
hood (colour-coded) and the 99% confidence level of the 
marginalised likelihood in the (6, log c) plane of param- 
eter space. The inclusion of large scale structure data 
and BAO data considerably tightens the constraints on 
features at small scales corresponding to values of b be- 
tween ~ 14.1 and ~ 14.4, while for larger values of b, i.e., 
features at larger scales, the contours remain roughly the 
same. 
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-4.5 ' ' ' ' 1 

14 14.2 14.4 14.6 14.8 15 
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FIG. 6: 99% confidence level contours for model B in the 
(b, logc) plane of parameter space with data sets 1 (top), 
2 (centre) and 3 (bottom). In these directions of parameter 
space, the likelihood function has a plateau towards vanishing 
step heights where the model reduces to the featureless m 2 cj> 2 
case, an excluded valley corresponding to large steps and a 
peak at 6 ~ 14.8. For the WMAP data alone we also find a 
second peak near b ~ 14.3. 



The feature at large scales (b ~ 14.8), on the 
other hand, remains untouched when we add the small 
scale data sets. For the generalised step model and 
data set 1, the large scale feature maximum likelihood 
point is at (6 = 14.8, c = 0.001, d = 0.02, Q b h 2 = 0.0216, 
fl c h 2 = 0.102, r = 0.11, n s = 0.952, logllO 10 ^] = 3.05 
and Hq = 72.7), which lies near the maximum of the 
marginalised ID likelihoods of the vanilla model in Fig. [5] 
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step), and, for T) > 1, not only £ 2 , but also higher order 
potential slow roll parameters will be non-negligible. 




FIG. 7: Marginalised (solid line) and mean likelihoods (dotted 
line) for parameter k s in the generalised step model, indicating 
at which wavelengths a feature is likely to happen. Top left: 
data set 1, top right: data set 2, bottom: data set 3. 



This is a further indication that the presence of a feature 
at large scales will not affect the estimates of the other 
parameters. 

Going from model B to the generalised step model will 
slightly improve the quality of the fits, yielding an extra 
A% 2 of 1-2. We did not expect a major improvement 
here, since the spectral tilt of the m 2 (j) 2 model lies fairly 
close to the best fit value of the vanilla model with a 
freely varying ng- 

In Fig. [Jj we show the marginalised and mean like- 
lihoods for the wavenumber k s of the perturbations 
that left the horizon when the inflaton field passed the 
step (i.e., at the moment when cf> — b). Again, we 
can see how the inclusion of the data sets sensitive to 
smaller scales reduces the evidence for a feature at scales 
> O(10~ 2 ) Mpc -1 . The difference between the mean and 
marginalised likelihoods is due to a volume effect: inte- 
gration over the low c plateau of the likelihood function 
tends to suppress peaks in the marginalised likelihood, 
which show up more clearly in the mean likelihood. 

Finally, we display constraints on the maximum val- 
ues of the slow roll parameters of the step function in 
Fig. [8] While the WMAP3 data alone is only sensitive to 
features up to a wavelength of ~ 10~ 2 Mpc -1 , the large 
scale structure data extends the sensitivity by almost a 
factor ten in k. We find fairly strong bounds on the maxi- 
mum value of e for the step function. In conjunction with 
Eq. ((55)) , this implies that the spectrum of tensor pertur- 
bations is unlikely to experience an oscillatory modula- 
tion like the scalar spectrum, since that would require e 
to be of order one. 

For the higher order slow roll parameters, values up to 
a few (for 77) and up to a few hundred (for £ 2 ) are still 
allowed. Note, however, that these bounds are param- 
eterisation dependent (they assume a tanh-form of the 




-2.5 -2 
Log 10 [k s /Mpc] 



FIG. 8: This plot shows the constraints on the peak values of 
the slow roll parameters during the step for the generalised 
step model. The thick line denotes the 99% confidence level 
for data set 3, the thin line corresponds to data set 1. 



VI. CONCLUSIONS 

We have analysed the dynamics of single-field infla- 
tion models with a step-like feature of small amplitude 
in the inflaton potential. Generically, the resulting spec- 
trum of scalar perturbations will resemble that of the 
stepless background model with a superimposed burst of 
oscillations whose shape is determined by the form of 
the step only. We have confronted the theoretical pre- 
dictions for the spectrum of a specific chaotic inflation 
model with a step with recent cosmological data to find 
out whether the data require the presence of such a fea- 
ture and whether it may actually bias the estimates of 
other cosmological parameters such as, e.g., the baryon 
density. We have also repeated the same analysis for a 
more empirical but less model dependent spectrum, such 
as might be expected from a step in an arbitrary infla- 
tionary background model. 

With a combination of different data sets, a large chunk 
of the step model parameter space can be ruled out, only 
spectra with a very modest oscillation amplitude are still 
consistent with observations. The BAO data, in partic- 
ular, prove to be a very sensitive probe for oscillating 
spectra. 

Compared to the 6 parameter "vanilla" cosmological 
model, using the most constraining data set, we find an 
improvement of the best fit \ 2 of about 5 for the chaotic 
inflation step model which comprises two extra parame- 
ters, and A\ 2 — 7 for the generalised step model, which 
has three extra parameters. 
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The vanilla model is a subset of the class of gener- 
alised step models for c — > 0. If c < 0(1O~ 5 ), the result- 
ing spectrum will be virtually indestiguishable from the 
vanilla spectrum. With our choice of priors, contours of 
greater than ~ 20% confidence level will contain parts of 
this vanilla region of parameter space. Hence, we cannot 
exclude the vanilla model at more than 20% confidence 
level. Reversing the argument, the present data do not 
show compelling evidence for requiring a spectrum with 
an oscillatory feature of the type discussed above. We 
expect that a more sophisticated model selection analy- 
sis along the lines of Refs. [H, IH) H3] would lead to a 
similar conclusion. 

The best fit region of parameter space consists of mod- 
els which show oscillations at wavelengths correspond- 
ing to multipoles t ~ O(10), where the temperature- 
temperature correlation data of the CMB shows some 
glitches. Interestingly, the time it would take the infla- 
ton field to traverse the step in these models is of the 
order of an e-folding, which is what one would expect for 
the time of a phase transition in more realistic multi-field 
models. 



Whether the glitches are just statistical flukes or stem 
from a physical effect, such as a feature in the inflaton 
potential, cannot be conclusively decided until we have 
better measurements of the E- and _B-mode polarisation 
spectra from experiments like PLANCK ,5l| or, in the 
more distant future, projects like the Inflation Probe [o2l ]. 
An additional consistency check can be provided by an 
analysis of the bispectrum of CMB fluctuations, since 
the interruption of slow-roll may also induce sizable non- 
Gaussianities [53[. 
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